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Abstract 

We predict the stabilities of a-graphynes and their boron nitride analogues(a-BNyne), which are 
considered as competitors of graphene and two-dimensional hexagonal BN. Based on first-principles 
plane wave method, we investigated the stability and structural transformations of these materials at 
different sizes using phonon dispersion calculations and ab-initio finite temperature, molecular dynam- 
ics simulations. Depending on the number of additional atoms in the edges between the corner atoms 
of the hexagons, n, both a-graphyne(?i) and a-BNyne(?i) are stable for even n, but unstable for odd 
n. a-graphyne(3) undergoes a structural transformation, where the symmetry of hexagons is broken. 
We present the structure optimized cohesive energies, electronic, magnetic and mechanical properties 
of stable structures. Our calculations reveal the existence of Dirac cones in the electronic structures 
of a-graphynes of all sizes, where the Fermi velocities decrease with increasing n. The electronic and 
magnetic properties of these structures are modified by hydrogenation. A single hydrogen vacancy 
renders a magnetic moment of one Bohr magneton. We finally present the properties of the bilayer a- 
graphyne and a-BNyne structures. We expect that these layered materials can function as frameworks 
in various chemical and electronic apphcations. 

Keywords: graphyne, hydrogenation, phonon modes, molecular dynamics, density functional the- 
ory 



Introduction 

The synthesis of graphene ^ has been an important 
turning point in the study of stable two dimen- 
sional (2D) monolayer structures. Graphene has 
drawn great attention due to its properties such as 
high chemical stability, mechanical strength and 
electric conductivity.^^ These unique electronic 
properties are mainly related to the Dirac cones 
present in its band structure where the upper and 
lower cones consisting of the conduction and the 
valence bands meet at a point on the Fermi level 
making graphene a semi-metal with zero band gap. 

Graphynes, similar to graphene, are two dimen- 
sional structures with the inclusion of single and 



triple bonded carbon atoms between the corner 
atoms of the honeycomb structure. Much earlier 
than the synthesis of single layer graphene, Baugh- 
man et al. ^ predicted various flakes or molecules 
of carbon atoms in the graphyne family as lay- 
ered phases using semi-empirical and empirical 
atom-atom potential calculations. These finite 
size nanostructures form from the combination of 
hexagons with other polygons containing sp^ and 
sp bonds. Based on first-principles plane wave 
calculation, Tongay et al^hme, predicted various 
stable ID, 2D and 3D periodic structures contain- 
ing carbon atomic chains. Among them, 2D pe- 
riodic a-graphyne have been revealed. Very re- 
cently, band structures of the graphyne/graphdyne 
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family with similar behaviors to that of single layer 
graphene were calculated,^ showing that neither 
the existence of hexagonal symmetry nor all atoms 
being chemically equivalent are prerequisites for 
the existence of Dirac point in the electronic struc- 
tures. 

Different types of the graphyne family are con- 
sidered as a new class of 2D materials in the fu- 
ture era of carbon allotropes7 Finite- size building 
blocks of these graphyne structures have already 
been synthesized which is an initial step towards 
extended structures .l^^^ Although extended (pe- 
riodic) 2D structures of Of-graphyne have not yet 
been synthesized, the synthesized finite size flakes 
building blocks are promising for future applica- 
tions. Theoretical studies and simulations have 
revealed that different members of the graphyne 
family can lead to interesting electronic applica- 
tions and they can be used to construct graphyne- 
based frameworks and nanotubes. It is also possi- 
ble to form the single layer, hexagonal boron ni- 
tride (h-BN) analogues of certain types of gra- 
phyne, which are called as BNyne throughout this 
manuscript.^ ^^'^^ Recently, effects of H, B and 
N doping into graphyne structures were also in- 
vestigated.!^ Additionally, the physical and chem- 
ical properties of graphyne and BNyne families 
and their future applications can be diversified by 
varying the sizes of rings and/or by functionalizing 
them with foreign ad-atoms. 

Although graphyne allotropes have been widely 
studied in previous works, their stabilities have 
been an open question. In this paper we inves- 
tigate the stabilities of 2D periodic oc-graphyne 
structures and their h-BN analogues by using first- 
principles calculations within density functional 
theory. An oc-graphyne structure can be obtained 
by placing n number of carbon atoms between the 
corner atoms of the hexagonal graphene structure. 
We calculate phonon modes and perform ab-initio, 
finite temperature molecular dynamic s(MD) sim- 
ulations to test the stabilities of different sized 
a-graphynes and a-BNynes. We investigate the 
electronic and mechanical properties of the stable 
structures. We find that the stabilities, the exis- 
tence of Dirac points at the Fermi level and the 
Fermi velocities of massless Fermions depend on 
the size of the ct-graphyne, such that n=even cases 
lead to stable structures with graphene like elec- 



tronic structures, whereas n=odd cases lead to in- 
stabilities and eventually undergo structural trans- 
formations. In addition, the mechanical strengths 
of a-graphynes of various sizes(n = 2 — 4) were 
explored and compared with those of other honey- 
comb like structures. We examine how the elec- 
tronic properties are affected by hydrogen atoms 
adsorbed at the corner carbon atoms of hexagons 
alternatingly from the top and bottom sides and 
thus formed a structure analogous to graphane.l^ 
In addition, we provide a comprehensive com- 
parison of cohesive energies with those of con- 
stituent ID and 2D allotropes of carbon. We fi- 
nally investigate the cohesive energies and elec- 
tronic properties of bilayer structures of these ma- 
terials. We also perform similar analysis for the h- 
BN analogues of graphynes, i.e. oc-BNynes. Our 
study revealed a new, non-hexagonal a-graphyne 
structure, where the symmetries associated with 
hexagons are broken. 

Method 

In our calculations we used the state-of-the-art 
first-principles plane-wave calculations within the 
density-functional theory combined with ab-initio, 
finite temperature MD calculations using projec- 
tor augmented wave potentials.!^ The numerical 
plane wave cal culati ons were performed using the 
VASP package.'^^'^ The exchange correlation po- 
tential was approximated by the generalized gra- 
dient a pprox imation with the van der Waals cor- 
rection.lSES] a plane-wave basis set with energy 
cutoff value of 600eV was used. The Brillouin 
zone (BZ) was sampled in the k-space within the 
Monkhorst-Pack scheme,'^ and the convergence 
of the total energies and magnetic moments with 
respect to the number of k-points were tested. In 
the self consistent total energy calculations and 
the MD simulations, the BZ was sampled by 
(17x17x1) mesh points in the k- space. The 
convergence criterion for energy was chosen as 
lO^^eV between two consecutive steps. In the ge- 
ometry relaxation and band structure calculations, 
the smearing value for all structures was taken as 
O.Oley. The pressure on each system was kept 
smaller than ~ IkBar per unit cell in the calcu- 
lations. In the ab-initio MD calculations, the time 
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Figure 1 : Atomic structure of a-Grapliyne and a-BNyne. (a) Scliematic representation of a-grapliyne(2) and tlie 
unit cell used to generate a-graphyne(«), where n is the number of carbon atoms placed between two carbon atoms 
located at the corners of the hexagon. Two corner atoms of the hexagon have a chain of n atoms between them, 
such that the unit cell contains 3n + 2 atoms, (b) Atomic structure of single layer, 2D a-graphyne(2). The dashed 
lines delineate the primitive unit cell. The optimized bond lengths are g\ = 1.39A and gl = 1.23A. The total charge 
density is shown within the unit cell. Carbon atoms are represented by brown balls, (c) Atomic structure of single 
layer, 2D a-BNyne(2) with blue and green balls representing N and B atoms, respectively. The optimized bond 
lengths are b\ = 1.42A, b2 = 1.25A and W = 1.44A. In the charge density plots, the isosurface value is taken as 
0.2 electron/ A^. 



step was taken as 2.5 fs and atomic velocities were 
renormalized to the temperature set at T = 500^ 
and T = 1000^ at every 40 time steps. In the 
MD stability tests, the simulations were run for 
5ps. The phonon dispersion curves were calcu- 
lated using plane wave methods as implemented 
in the PWSCF package.l^^l 

a-Graphyne and a-BNyne 

Structure 

We begin our analysis by investigating the struc- 
tural and geometrical properties of ct-graphynes 
at different {n) sizes. As shown in [TJa), a- 
graphyne(n) structure has a hexagonal unit cell 
like graphene, but there are 3n + 2 carbon atoms 
in its unitcell instead of 2. In other words, we have 
the normal graphene-like structure with chains'^ 
consisting of n carbon atoms between the corner 
atoms of the hexagon on each side. Contrary to 
graphene, all of the C-C bonds are not equiva- 
lent to each other in these structures, but there ex- 
ists bonds with different lengths and charge den- 
sities as shown in[T]^b). Variances of bonds origi- 
nate from different types of bonding, sp^ and sp^ 
i+Tl) between carbon atoms. Therefore, the carbon 
atoms are no longer chemically equivalent to each 
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Figure 2: Calculated phonon bands, (a) Graphene. (b) 
a-Graphyne with n =2 and n =4. (c) Single layer h- 
BN. (d) a-BNyne with n = 2 and n =4. The disper- 
sion curves for n=2 have totally positive phonon modes 
which is an indication of their stability. On the other 
hand, «=4 cases have modes with imaginary frequen- 
cies, which are marked with the shaded regions and 
will be discussed in the text. Phonon bands of unsta- 
ble structures, such as « = 1 and « = 3 are not shown. 
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a-Graphyne[n) 




Figure 3: Snapshots of the MD simulations performed for 5ps a\T = \000K. (a) a-graphyne(n) n =1-4 structures. 
The structures are stabie for n =2 and n =4, although buckied in the vertical plane. On the other hand, n = I case 
breaks into carbon atomic strings, and hence are totally unstable. a-graphyne(3) undergoes a structural transforma- 
tion, whereby it acquires stability by changing the number of C atoms to « = 2, n = 3 and n = 4 \n the adjacent edges 
of hexagon, (b) Sections of a-BNyne(?i). Both n =2 and n =4 cases remain stable during MD simulations, n =1 
and n =3 cases are missed, since a-BNyne(n) cannot be formed with odd n. Note that only a single ring from the 
periodic structures are shown here. 
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Figure 4: (a) Top and side views of the optimized 
geometry of the structurally transformed, stable a- 
graphyne(3). The planar, unstable a-graphyne(3) un- 
dergoes a structural transformation and acquires stabil- 
ity by breaking the symmetry of hexagons and eventu- 
ally by buckling. The adjacent edges of the hexagon has 
n=2, 3 and 4 carbon atoms, respectively, (b) Electronic 
structures of the buckled a-graphyne(3). The buckling 
leads to a gap opening of 0.42eV. 



other as they were in the case of graphene. 

The h-BN analogue of this structure, namely 
a-BNyne, can be easily obtained by replacing 
the carbon atoms with B and N consecutively. 
The electronegativity numbers of B and N atoms 
are 2.0 and 3.0, respectively according to Paul- 
ing's electronegativity scale. Hence the elec- 
tronic charge is transferred from B to adjacent N 
in the ionic bonds of a-BNyne. Here, there are 
three different types of bonds; two at the corners 
of the hexagon, one at the center of edges. For B- 
N bonds at the corners, one is the case where N is 
at the corner of the hexagon and B in the edge as 
denoted hy bl in[TJc), and the other is the opposite 
of this and denoted by b3. Note that, if we want to 
preserve the 2D-BN hexagonal structure such that 
the corner atoms of the hexagons are consecutively 
B and N, we must have even number of atoms be- 
tween the comer atoms. That is, we can only have 
n =even scenarios for the a-BNyne(n) cases. Odd 
n values would either lead to same type of atoms 
on the corners of the hexagons or two identical 
atoms(B-B or N-N) next to each other. For this 
reason, we restrict ourselves with even n values 
for a-BNyne(n). The periodic geometries of these 
structures, and the unit cells used to generate them 
are shown in[TJb) and[T]^c) for a-graphyne(2) and 
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a-BNyne(2) cases, respectively. 

Stability 

Up to now, finite building blocks of graphyne 
allotropes'^'^ and also the segments of atomic 
chains consisting of n carbon atomsl^^^ have 
been synthesized, which may be taken as an in- 
dication of the stabilities of their extended struc- 
tures. Here, the main issue is whether 2D single 
layer, periodic structure of oc-graphyne is stable 
or not. As an initial step we calculated the co- 
hesive energies of these structures, which is de- 
fined as the energy required to form separated neu- 
tral atoms in their ground electronic state from the 
condensed state at 0^ and latm!^ Cohesive en- 
ergies per atom of the carbon allotropes (namely 
graphene, a-graphyne(/i) and cumulene or infinite 
carbon atomic chain ^'^) are obtained using the ex- 
pression 

Ec = {pEj [C] - Et [Callotrope]}/p, (D 

in terms of the optimized total energies of single 
carbon atom in its magnetic ground state Et[C\, 
and of the related carbon allotrope EjlCaiiotrope], 
with p being the number of carbon atoms in the 
unit cell. In the case of BN structures, the cohesive 
energy is calculated for per B-N pair using 

Ec = {qET[B\ +qET[N] -ET[BNaiiotrope]} / q (2) 

in terms of optimized total energies of single 
B atom Et[B\, N atom Et[N] and of the related 
BN allotrope, with q being the number of B-N 
pairs in the unit cell. Here (Q^)aUotrope stands 
for one of the structures, such as BN chain, h- 
BN and a-BNynes. The calculated cohesive en- 
ergies of graphene, a-graphyne(n), cumulene and 
their h-BN analogues are given in Table 1. We 
note that the cohesive energies of a-graphyne(n) 
and a-BNyne(n) are smaller than their parent, sin- 
gle layer graphene and h-BN honeycomb struc- 
tures, respectively. Also expectantly, the cohe- 
sive energies of a-graphyne(n) and -BNyne(n) de- 
crease with increasing n, since the cohesive en- 
ergy of the linear atomic chain (either cumulene 
or BN) is smaller than the corresponding 2D sin- 
gle layer, honeycomb structure (respectively, ei- 



ther graphene or h-BN). Like graphene and carbon 
atomic chains, the calculated cohesive energies 
suggest that a-graphyne structures correspond to 
local minima on the Born-Oppenheimer surface. 
Accordingly, even the network of carbon atomic 
chains consisting of diverse number of atoms n, 
which are connected by three folded 5/7^-orbitals 
at the nodes can be in local minima.'^ In principle, 
these networks can form either multiply connected 
polygons(rings) or tree-like structure. Here we 
note that the cohesive energies of a-graphyne(n) 
for n =1 and n =3 calculated by restricting the 
structure in a plane are found smaller than that of 
n =2 and n =4. As clarified in the following anal- 
ysis of stability using phonon calculations, n =odd 
number structures are unstable despite their posi- 
tive cohesive energy. Here, smaller Ec for n =1 
and n =3 relative to that of n =2 and n =4 in- 
dicate instability, but other analysis like calcula- 
tion of phonon dispersion and MD simulations are 
needed for a reliable test of stability. 

Table 1 : Cohesive energies(£'c) of optimized grapliene, 
a-grapliyne(«) n=^-4, cumulene structures and their h- 
BN analogues. The cohesive energies of the BN struc- 
tures are given per number of B-N pairs. 



Structure 


Ec (eV) 


Graphene 


8.11 


Carbon chain 


7.10 


a-graphyne(l) 


6.21 


a-graphyne(2) 


7.35 


a-graphyne(3) 


6.82 


a-graphyne(4) 


7.12 


2D BN 


14.52 


BN chain 


12.80 


a-BNyne(2) 


13.26 


a-BNyne(4) 


12.90 



A state-of the art calculations of phonon fre- 
quencies of oc-graphynes and a-BNynes are car- 
ried out for all modes as a function of the k-points 
in the BZ. When the calculated frequencies of 
all the modes are positive, the structure is iden- 
tified to be stable. We note that the constituent 
allotropes, namely graphene, h-BN, carbon and 
BN chains, have positive frequencies for all modes 
in their BZ, and hence they are stable. As seen 
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in |2} all calculated frequencies of 0C-graphyne(2) 
and a-BNyne(2) are positive. Moreover, due to 
the segments of atomic carbon and BN chains in 
a-graphyne(2) and a-BNyne(2), both structures 
have phonon modes with frequencies higher than 
those of single layer graphene and h-BN. As a mat- 
ter of fact, the maximum frequency of the longitu- 
dinal optic mode (LO) of carbon chain'^can be as 
high as ~ 2400 cm^^ 

The situation is seemingly different for a- 
graphyne(4) and a-BNyne(4), since some of 
the acoustic modes have imaginary frequencies. 
These modes shaded out in [2fb) and (d) cor- 
respond to imaginary frequencies and normally 
would indicate instabilities. However, the imag- 
inary frequencies may arise also as an artifact of 
numerical calculations. To obtain the real frequen- 
cies of soft modes in an open structure comprising 
many atoms, one need to perform calculations by 
taking into account the distant neighbors and use 
very high numerical accuracies. Thus, we be- 
lieve that the calculated imaginary frequencies of 
soft modes in this case are artifacts of numerical 
phonon calculations since the segments of chains 
consisting of even carbon atoms have been found 
to be stablJ^ and they become even more stabi- 
lized if they are connected to carbon atoms with 
three folded sp^ orbitals. Also, it should be noted 
that all a-graphyne(n) and a-BNyne(n) with large 
n can be stabilized at finite-size, since long wave- 
length acoustical modes are discarded. 

In order to further investigate the issue of stabil- 
ity, we performed MD simulations at T = 1000^ 
for 5ps. The atomic structures obtained after the 
MD calculations are presented in [3| It turns out 
that for n =2 and n =4, a-graphyne(n) and a- 
BNyne(n) remain stable after 5ps of MD simula- 
tion even at high temperatures. Even though 5ps 
is a short time interval, it is long enough for ab- 
initio calculations at the temperature as high as 
T = 1000^ to provide evidence for stability. The 
analysis of atomic structures also suggests that a- 
graphyne(4) and a-BNyne(4) prefer to buckle in 
the vertical direction. Also in some cases, when 
specific structures are unstable in planar geometry 
they might undergo structural transformations to 
attain stability. Silicene {i.e. Si in graphene struc- 
ture) is a crucial example for stability gained by 
buckling. -^^ Here we assure that the buckling oc- 



curs due high temperature for the following rea- 
sons: (i) The optimization of buckled structure 
at T=OK ended with planar structure; (ii) similar 
buckling effects were also observed at segments of 
carbon atomic chains self-assembled on graphene 
and h-BN which bowed at high temperatures.'^^^®' 
It was found that a-graphyne(l) is totally unsta- 
ble and is dissociated into carbon atomic strings. 
a-graphyne(3) presents an interesting situation; it 
undergoes a structural transformation and acquires 
stability by changing the number of carbon atoms 
to n =2, n =3 and n =4 in the adjacent edges of 
hexagon as shown in|4j This structural transforma- 
tion is derived to maintain the proper bond order of 
finite size carbon chains constituting the edges of 
hexagon.f^ While the carbon atoms at the corners 
are always forming single bonds with the adjacent 
carbon atoms, the second atom from the corner 
by itself has to make a triple bond with the ad- 
jacent carbon atom, which is at the other side in 
the edge of the hexagon. At the end, the correct 
bond order of carbon atoms are preserved. How- 
ever, this structural transformation, which modi- 
fies the geometric structure and breaks the symme- 
try of hexagons, reveals a new buckled alio trope 
carbon atom in 2D with a band gap of 0.42eV. 




Figure 5: (a)a— Graphyne(2) structure in rectangular 
unit cell with its lattice constants Qv and Sx and Sy 
are the strains in x and y directions, respectively, (b) 3D 
plot of the energy values corresponding to different ax 
and ay values. 



Mechanical properties 

Having found the stable a-graphyne and a-BNyne 
structures, we next calculate their mechanical 
strengths. A common way of expressing the me- 
chanical properties of two dimensional materials is 
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Figure 6: Electronic band structures of a-grapliyne(?i) for « = 1,2,3 and 4. All of the band structures contain Dirac 
points, while they are shifted above the Fermi level for « = 1 and 3. « = 1 and 3 cases also have Dirac points away 
from the high symmetry ^— point. The zero of energy is set to the Fermi level. Note that, the electronic structures of 
the n = 1 and n = 3 structures are presented for the sake of completeness although they are unstable in the planar 
configuration, but « = 3 structure acquires stability in the buckled geometry as discussed in the text. 
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Figure 7: Electronic band structures of stable a- 
BNyne(?i) for « = 2 and 4. Note that as n increases, 
the band gap decreases. The maximum energy of the 
valence band is set to zero. 



to calculate their in-plane stiffness, Poisson's ratio 
and Young's modulus values. For this purpose, we 
used a rectangular unit cell in the xy-plane and ap- 
plied tension in both x and y directions, as shown 
in[5]^a). We varied the lattice constants, and 
between ±0.03% of the optimized values and cal- 
culated the energy value for each grid point ob- 
tained. In the end, we have obtained energy values 
for 225 grid points, which are plotted in [Sj^b) for 
a-graphyne(2). The strain energy E^, at each point 
is calculated by subtracting the total energy at that 
point from the equilibrium total energy. It has been 
shown that this method of calculating the strain 
energy provides reliable predictions for graphene, 
and 6-6-12 graphyne.S^flsj 

The in-plane stiffness, which is commonly used 



measure of strength for 2D materials can be ex- 
pressed as 



C 



1 d'^E, 



(3) 



where Es is the strain energy, Aq is the equilib- 
rium area and e is the axial strain calculated by 
l^a^.y! a^y, a being the lattice constant in the xoiy 
direction. 

The in-plane stiffness values of a-graphyne(2), 
a-graphyne(4), a-BNyne(2) and a-BNyne(4) 
were calculated as 21N/m, l6N/m, l9N/m, and 
IAN /m, respectively. These values are much lower 
than the in- plane stiffness of graphene, which is 
~ 340A'^/m.S^'^ As seen from these results, the 
implementation of new atoms between atoms at 
the corners of hexagons decreases the mechanical 
strength of graphene dramatically, which is a di- 
rect consequence of the decrease in the average co- 
ordination number of carbon atoms. The Poisson's 
ratio, which is defined as the ratio of the trans- 
verse strain to the axial strain, v = —£trans/£axiah 
was calculated as 0.88, 0.86, 0.89 and 0.85 for 
a-graphyne(2), a-graphyne(4), a-BNyne(2) and 
a-BNyne(4), respectively. By assuming an equiv- 
alent thickness with graphene, the Young's modu- 
lus values were calculated, respectively as 6lGPa, 
ASGPa, 52GPa and A2GPa for these structures. 
We also note that in-plane stiffness, Poisson's ratio 
and Young's modulus values of a-graphyne(n) and 
a-BNyne(/?) decrease with increasing n. Also, the 
values calculated for a-BNyne(n) are lower than 
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those calculated for a-graphyne for each case with 
equivalent n. 

Electronic structure 

Earlier, Tongay et al. ^ showed that the electronic 
structure of a-graphyne(2) with two bands cross- 
ing the Fermi level at K- and ^'-points of the BZ 
is similar to that of graphene. The presence of 
Dirac cones in the band structure of ct-graphynes, 
as well as /3— graphyne and 6-6-12 graphyne have 
also been recently reported,^-* and it was demon- 
strated that the existence of Dirac points and cones 
is not a unique property of graphene. The crossing 
of bands at the Fermi level and the formation of 
Dirac cones were also investigated for other struc- 
tures.*^ It was pointed out that bare structures, 
large defects and ad-atoms on graphene can have 
Dirac cones if their periodic patterns comply with 
a specific symmetry. Here, we present the elec- 
tronic energy band structure of a— graphyne(n) for 
n =1,2,3 and 4 in[6j The existence of Dirac points 
is also seen here for all of these graphynes. Note 
that, the Dirac point lies at the Fermi level for the 
stable n=2 and n=4 structures, whereas it is shifted 
above from the Fermi level for the unstable n=\ 
and n=3. We note that high electron density at 
the Fermi level of a— graphyne(n) for n=\ and n=3 
can be attributed to their instability. The electronic 
energy near the K-point of the BZ is linear with 
respect to q=K-k, which leads 

E{±) = ±hvfq + 0[{q/k)\ (4) 

where v^- is the Fermi velocity and K is the wave 
vector corresponding to K- and ^'-points of the 
BZ. Then, for the stable a-graphyne(2) and a- 
graphyne(4) structures, the first derivatives of their 
n bands near the ^-point of the BZ were calcu- 
lated as 29.4 eVA and 23.4 eVA , which have the 
same order of magnitude obtained for graphene, 
34.6 eVA . We also estimate the Fermi velocities 
as vp ~ 8.3 X 10^ m/s for graphene, v/r ~ 7.1 x 10^ 
m/s for a-graphyne(2) and vf ~ 5.6 x 10^ m/s for 
a-graphyne(4). Accordingly, the Fermi veloci- 
ties at K- and K -points decrease with increasing 
n. Noting the electron-hole symmetry of Eq(4) 
near ^-point, a-graphyne(2) and a-graphyne(4) 
are found to be ambipolar. As discussed in Sec. 



Ill B, a-graphyne(3) having metallic state undergo 
structural transformation where hexagons forming 
a honeycomb structure transforms to rectangle like 
structures as shown in the third panel of |3j When 
planar, this structure attains high density of states 
at the Fermi level. However, through the buck- 
ling of the atomic structure it undergoes a metal- 
insulator transition by opening an indirect band 
gap of 0.42 eV (See ^. 

On the other hand, reminiscent of the electronic 
structure of 2D single layer h-BN, the band struc- 
tures of a-BNyne(2) and a-BNyne(4) have wide 
band gaps as shown in |7j It can be seen that the 
band gap decreases with increasing values of n. 
The lowest conduction band and the highest va- 
lence band states originate from B-p^ and N-p^ or- 
bitals. 

Hydrogenation 

As an immediate application a-graphyne(n) and 
a-BNyne(n), one may consider their chemical 
conversion through the coverage of H (hydrogena- 
tion), F (fluorination), CI (chlorination). Here we 
present our results regarding the hydrogenation of 
a-graphyne(2) and a-BNyne(2). The hydrogena- 
tion of graphene which produces graphane, and 
the dehydrogeneation of graphane are already well 
known process. '^^'^ Subsequently, hydrogenation 
of single layer BN was also theoretically studied 
which lead into to BN-phane.^^ In a similar way, 
we attach single hydrogen atoms to the corner car- 
bon atoms of a-graphyne(2) and B and N atoms 
at the alternating comers of a-BNyne(2) from top 
and bottom alternatingly, as illustrated in [Sj^a-b). 
Upon hydrogenation, both a-graphyne(2) and a- 
BNyne(2) structures relax into a buckled geome- 
try with buckling distances of 0.92 A and 1.03 A. 
These buckled geometries are similar to the buck- 
led geometries of graphane and BN-phane. How- 
ever, the buckling increases from 0.45 A to 0.92 
A as we go from graphane to hydrogenated a- 
graphyne(2), as a result of the increased length of 
the edges of hexagons. Hydrogenation also alters 
the electronic structures of these systems, as pre- 
sented in [8]^c-d). Consequently, the semimetalic 
a-graphyne(2) structure attains a wide band gap of 
5.2 eV upon hydrogenation. While hydrogenated 
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a-graphyne(2) has nonmagnetic ground state, a 
single hydrogen vacancy renders a magnetic mo- 
ment of 1 jlB 

Hydrogenated a-BNyne(2) also has a wide band 
gap of 3.3eV. The band gap opening effect of 
hydrogenation on monolayer structures is in ac- 
cordance with prior calculations on graphanel^ 
and hydrogenated BN,^^ which have band gaps of 
3.5eV and 4.8eV, respectively. 




TMKF FMKr 

Figure 8: (a-b) Hydrogenated a-graphyne(2) and a- 
BNyne(2). Top and side views of tine optimized liy- 
drogenated structures are sliown by tine ball and stick 
models, where C, N, B and H atoms are represented 
by brown, blue, green and red balls, respectively, (c- 
d) Electronic structures of hydrogenated a-graphyne(2) 
and a-BNyne(2). Note that, in contrast to the semimet- 
alic a-graphyne(2), the hydrogenated structure has a 
wide band gap of 5.2eV. The maximum energy of the 
valence band is set to zero. 



Bilayer Structures 

Here we address the question of whether a- 
graphyne and a-BNyne can form layered struc- 
tures similar to graphite and hexagonal BN, or 
not. We place their bilayers as shown in fig- 



ure |9] and explore equilibrium geometries, bind- 
ing energies and electronic structures. We begin 
the analysis by placing two single layers of a- 
graphyne(2) sheets on top of each other with AA 
stacking {i.e. hexagons in both layers face each 
other) and AB stacking (i.e. first layer is shifted 
laterally to the centers of hexagons in the sec- 
ond layer) geometries. We alter the interlay er dis- 
tances until we achieve the minimum energy val- 
ues. The calculated minimum energy geometries 
indicate that AB type of stacking is more favor- 
able than AA type by 46meV. This is a behavior 
similar to bilayer graphene structure or graphite. 
The optimized interlayer distance is calculated as 
3.12A, which is less than the interlayer distance 
of graphite, 3. 35 A. This is mainly related to the 
less dense arrangement of carbon atoms on the 
graphyne surface as compared to the graphene, 
which results in lower surface energy and hence 
closer equilibrium distance. Similar results were 
also found in a previous interlayer distance and 
stacking studies made on different graphyne al- 
lotropes.^^ 

A convenient procedure for calculating the inter- 
layer binding energy for layered structures is sub- 
tracting the minimum energy of the bilayer struc- 
ture from the sum of energies of separated indi- 
vidual layers."^ Applying this method, we calcu- 
late the binding energy of bilayer a-graphyne(2) 
as 220meV (21 .5meV / atom) . We repeat the same 
analysis for bilayer a-BNyne(2), which also has 
an AB type double layer geometry as shown in 
the second column of |9]^a). For this case, the en- 
ergy difference between AB stacking and the AA 
stacking is in favor of AB by lOmeV, the inter- 
layer distance is 2. 9 A and the binding energy is 
l2SmeV{l6meV /atom). The calculated binding 
energies of bilayer a-graphyne(2) and bilayer a- 
BNyne(2) are in the range of the calculated bind- 
ing energies of bilayer graphene (ASmeV /atom) 
and bilayer h-BN (31. 5meV /atom), respectively. 
The variations of total energies around the min- 
imum energy values as a function of the inter- 
layer distances are shown in ^h). Like graphite 
and layered BN, small interlayer binding ener- 
gies consist of mainly from the van der Waals in- 
teraction with even smaller chemical interaction 
component. Nevertheless, such a weak binding 
is enough to maintain the bilayer structure at low 
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temperature. Additionally, 3D layered structures 
of Ct-graphyne and a-BNyne can form at room 
temperature. These arguments are corroborated 
by the MD simulations performed at different tem- 
peratures showing that the bilayer structures form 
stable interlayer binding near room temperature, 
but the layers move away from each other at high 
temperatures {T = 500K and 1000^). These re- 
sults imply that multilayers of a-graphyne and a- 
BNyne, even their 3D layered structures can form. 



Biliyer (i-Graphyne{3) 
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Figure 9: Bilayer a-grapliyne(2) and its BN analogue 
bilayer a-BNyne(2) are shown in columns / and ii, re- 
spectively, (a) Top view of the optimized two layer struc- 
tures. Both bilayer a-graphyne and a-BNyne have AB 
type of stacking geometry, which is more favorable than 
the AA stacking. In the ball and stick model C, B and 
N atoms, respectively are represented by brown, green 
and blue balls and the all of the atoms in the bottom 
layer are shown in gray balls, (b) Variation of energy 
as a function of the layer-layer distance, (c) Electronic 
band structures of a-graphyne(2) and a-BNyne(2). 

We finally calculate the electronic structures of 
these double layered structures. The band struc- 
ture of bilayer a-graphyne(2) is shown in ^c). 
Note that, the bands are split due to the cou- 
plings between layers and the numbers of energy 



bands are doubled as compared to the single a- 
graphyne(2) sheet. Notably, the bands are no 
longer linear around the ^-point, but parabolic. 
This kind of structure is reminiscent of the elec- 
tronic structure of bilayer graphene or graphite.!^ 
As opposed to single layer, the bands no longer 
touch at the K point but there is a small band gap 
of lOmeV. The electronic structure for the bilayer 
a-BNyne(2) is also presented in|9jc). Similar to 
the single layer, the double layered structure also 
has a wide band gap. 

Conclusion 

In conclusion, we investigated the cohesion, struc- 
tural stabilities, electronic structures and mechan- 
ical properties and functionalization by adatoms 
of a-graphyne(/?) structures along with their h- 
BN analogues a-BNyne(R). By performing both 
phonon frequency and finite temperature molec- 
ular dynamics analysis for n = 1 — 4 cases, we 
showed that both a-graphyne and a-BNyne are 
stable structures for even n, but unstable for odd 
n. Interestingly, a-graphyne(3) undergoes a struc- 
tural transformation to acquire stability, whereby 
the symmetry of hexagons forming honeycomb 
structure is broken. Thus, our study clarified the 
question whether 2D periodic a-graphyne and a- 
BNyne structure can be stable or not. We also cal- 
culated the electronic structures for each of these 
materials and showed that a-graphyne(n) struc- 
tures having Dirac cones are ambipolar and their 
Fermi velocities decrease with increasing n. It 
is also implied that, all atoms being chemically 
equivalent is not a prerequisite for the existence 
of Dirac cones in the electronic structure. Upon 
hydrogenation, the Dirac cones are replaced by a 
large band gap. Since the formation of a single 
hydrogen vacancy renders a magnetic moment of 
1 /is, magnetic nanomaterials can be designed by 
creation of domains of hydrogen vacancies. Our 
calculations of mechanical properties revealed a- 
graphyne and a-BNyne are not as stiff as graphene 
and the single layer h-BN, but they are strong 
enough to sustain the technological applications. 

We finally showed that it is also possible to have 
double layers of a-graphyne and a-BNyne struc- 
tures. Both of these bilayer structures have AB 
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type of stacking. It was found that the electronic 
structure of bilayer oc-graphyne has a gap opening 
of lOmeV at the ^— point of the BZ as opposed to 
single layer a-graphyne(2). On the other hand, the 
electronic structure of bilayer a-BNyne(2) con- 
tains a wide band gap, similar to single layer BN 
structure. 

Briefly, we demonstrated that a-graphyne and 
a-BNyne nanostructures can form stable and 
durable 2D extended structures with interesting 
chemical and physical properties, which are scaled 
by n. We believe that, these stable 2D carbon and 
BN allotropes will attract interest because of their 
unique properties in the near future. In particular, 
they can be utilized as structural frameworks for 
various spintronic^^ and chemical applications. 
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